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Abstract 



We use two renormalization techniques, Effective Field Theory and the Similarity Renor- 
malization Group, to solve simple Schrodinger equations with delta-function potentials in one 
and two dimensions. The familiar one-dimensional delta-function does not require renormal- 
ization, but it provides the simplest example of a local interaction that can be replaced by 
a sequence of effective cutoff interactions that produce controllable power-law errors. The 
two-dimensional delta- function leads to logarithmic divergences, dimensional transmutation 
and asymptotic freedom, providing an example of some of the most important renormaliza- 
tion problems in gauge field theories. We concentrate on the power-law analysis of errors in 
low-energy observables. The power-law suppression of the effects of irrelevant operators is 
critical to the success of field theory, and understanding them turns renormalization group 
techniques into powerful predictive tools for complicated problems where exact solutions are 
not available. 



1 Introduction 

The theory of renormalization is one of the greatest successes of physics in the twentieth 
century Q. Early attempts to marry quantum mechanics with special relativity naturally 
led to the use of local interactions, maintaining the causal paradigm that has driven physics 
for centuries while avoiding signals that propagate at speeds greater than that of light. 
Early encounters with perturbative divergences led some of the best theorists in the world 
to question the foundations of quantum mechanics, and eventual successes at fitting precise 
atomic experimental data led to the universal acceptance of renormalization recipes that 
were acknowledged to make little sense [g|. Initially the perturbative renormalization of 
QED required theorists to match perturbative expansions in powers of a bare and physical 
electronic charge ||, but the bare charge clearly diverges logarithmically in QED and the 
success of an expansion in powers of such a coupling was mysterious at best || . 

The first steps towards making sense of renormalization theory were taken in the 1950's 
with the invention of the perturbative renormalization group || ^, ^, ||, [| , although serious 
investigators found the theory was still plagued by non-convergent sums because QED is not 
asymptotically free. The development of Wilson's renormalization group formalism |l^, |1J 
and the discovery of asymptotic freedom [|l^] allowed physicists to produce a logically rea- 
sonable picture of renormalization in which perturbative expansions at any high energy scale 
can be matched with one another, with no necessity to deal with intermediate expansions 
in powers of a large parameter. But in all theories of interest the running coupling becomes 
large at some scale and the lattice [|l3| provides a tool for attacking gauge theories non- 
perturbatively. This shifted the problems to achieving a continuum limit in which the lattice 
spacing can be taken to zero implicitly and set the stage for the development of a standard 
model that rests heavily on asymptotically free field theories. 

Although this story took a dramatic turn with the invention of string theories, its main 
themes survive in essentially all predictive field theory calculations and the control the renor- 
malization group provides has made it possible for physicists to accept the standard model as 
an effective theory that is accurate over a very large range of scales that covers all currently 
experimentally accessible energies. It has also led to the acceptance of the development of 
effective field theories as a fundamental endeavor that goes far beyond the apparently arbi- 
trary introduction of non-local potentials, which has persistently undermined broad interest 
in attempts to model complicated interactions such as the strong nuclear interaction. 

We make no attempt to survey this history, nor are we able to provide a complete review of 
the two methods illustrated in this article. We take advantage of the fact that the divergences 
that provided the original impetus leading to renormalization theory result entirely from the 
use of local interactions. To understand the most important aspects of renormalization theory 
requires only a background in nonrelativistic quantum mechanics, because as has been long 
known the divergences of field theory are directly encountered when one tries to impose 
locality on the Schrodinger equation. In this case the only available interactions that are 
local at all scales are delta functions and derivatives of delta functions. These divergences 
can be regulated by the introduction of a cutoff, and the artificial effects of this cutoff 
must be removed by renormalization. The simplicity of the one-body Schrodinger equation 
makes it possible to renormalize the theory exactly, disentangling the effects of locality from 
the complicated many-body effects and symmetry constraints encountered in realistic field 
theories. There is a large literature on the subject [14]-[31], largely pedagogical. 



Effective Field Theory (EFT) is primarily used to replace complicated 'fundamental' theo- 



ries such as QCD with simpler theories such as chiral perturbation theory [32 1 . These simpler 
effective theories are to be used only at 'low' energy scales at which the fundamental degrees 
of freedom and interactions cannot be resolved, and the scaling analysis that we illustrate 



makes this limitation clear. EFT |33|- |39[| is motivated by Wilson's work on the renormaliza- 
tion group jn| [llj; however, practitioners typically employ only basic renormalization group 
tools. 

EFT and essentially all renormalization techniques rely on the assumption that physics 
at a low energy scale and/or long distance scale is insensitive to the details of the underlying 
physics at a high energy scale and/or short distance scale. Microscopic degrees of freedom 
and their interactions can be replaced by effective macroscopic degrees of freedom and their 
effective interactions. Moreover, these effective degrees of freedom can be treated as if they 
were point-like and their effective interactions are local at the scale of interest. These as- 
sumptions have clearly been critical through the entire evolution of physics. However, the 
detailed procedure for producing effective interactions and the understanding of how errors 
depend on the parameters in these interactions are relatively new. 

One of the most interesting recent applications of EFT is to nuclear physics and the 



two-nucleon problem in particular |40|- |75jl . Traditionally nuclear theorists have developed 
nonlocal two-body interactions, sometimes motivated by meson exchange models, that are 
adjusted to fit scattering phase shifts and properties of the deuteron. These nonlocal inter- 
actions are typically parameterized by a few masses and couplings, but there has been no 
systematic, logically simple procedure for improving interactions and many nonlocal interac- 
tions work equally well. In EFT the cutoff is regarded as artificial and parameters are allowed 
to run with the cutoff to remove its effects as much as possible. There is a systematic pro- 
cedure for improving the effective interactions by adding operators, and the goal of relating 
the complicated nucleon-nucleon interaction to QCD may become much more realistic. 



Recent work to extend this approach to deal with the three- nucleon problem [76, 77] and 
to explicitly include mesons has uncovered problems that go beyond a straightforward appli- 
cation of EFT. The fact that the two-nucleon interaction is an irrelevant operator, but the 
scattering phases shifts are large, has led to an understanding that naive dimensional scaling 
is violated. There is a non-free fixed point that contains this interaction and presumably 



alters the scaling properties of all operators [47, 54, [57], 59]. Condensed matter theorists 
have long appreciated the importance of understanding the scaling behavior of interactions 
around non-Gaussian fixed points, but there are few field theory examples that are relevant 
to nuclear and particle physics. The study of effective nuclear field theories may offer such 
examples and then drive interesting new developments in EFT. 

The similarity renormalization group (SRG) is a very recent development invented by 



Stan Glazek and Ken Wilson [^, 79], and independently by Franz Wegner ]pQ ]. Most of 



the applications of the SRG are in light-front field theory [81[-[|95|], but there are a growing 



number of applications in other areas |96]-| 108 |. 



In EFT a cutoff is introduced along with a fixed number of effective interactions. The 
resultant hamiltonian is used to compute low energy observables, and the strengths (i.e., 
couplings) of the effective interactions are adjusted to fit a subset of the low energy data. 
It is sometimes possible to explicitly derive equations that govern how these couplings 
evolve as the cutoff changes, an approach that parallels the Gell-Mann-Low approach to 
the renormalization group [0], and by studying these equations the scaling behavior of these 
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Figure 1: Two ways to run a cutoff on free energy. In (a) a cutoff on the magnitude of the energy 
is lowered from the solid to the dashed lines, with problems resulting from the removed shaded 
region. In (b) a cutoff on how far off diagonal matrix elements appear is lowered from the dashed 
to the solid lines. 



couplings can be understood. In the SRG, as in Wilson's original renormalization group 
formalism [106, |HCf| , transformations that explicitly run the cutoff are developed. These 
transformations are the group elements that give the renormalization group its name. 



In his earliest work [109, 110] Wilson exploited a transformation originally invented by 



Claude Bloch [111]. It uses a cutoff on the states themselves, and as the cutoff is lowered, 
states are removed from the Hilbert space. If the hamiltonian is viewed as a matrix, these 
cutoffs can be seen as limiting the size of this matrix and the transformation reduces this 
size, as illustrated in Fig. la. Wilson introduced a rescaling operation to allow transformed 
hamiltonians to be compared with initial hamiltonians, despite the fact that they act in 
different spaces; however, the Bloch transformation is ill-defined and even in perturbation 
theory it leads to artificial divergences. These divergences come from the small energy dif- 
ferences between states retained and states removed by the transformation, and they appear 
in the form of small energy denominators in the perturbative expansion of the transformed 
hamiltonian. These small energy denominator problems led Wilson to abandon the hamilto- 
nian formulation of field theory in favor of path integral formulations, but the virtues of the 
hamiltonian formulation over the path integral formulation for many problems remains. 

The breakthrough provided by the SRG is that the transformations are typically unitary, 
making them well-defined, and they run a cutoff on energy differences rather than on indi- 
vidual states, as illustrated in Fig. lb. Again viewing the hamiltonian as a large matrix, 
these cutoffs limit the off-diagonal matrix elements and as they are reduced the hamiltonian 



is forced towards diagonal form. The perturbative expansion for transformed hamiltonians 
contains no small energy denominators, so the expansion breaks down only when interactions 
become sufficiently strong, in which case perturbation theory should fail in any case. 

Although the SRG has not yet been applied to a wide range of problems, it may be an 
important new tool both for attacking field theories and non-relativistic many-body problems. 
As we will illustrate, the effective interactions that arise in the SRG are nearly identical to 
those encountered in EFT; however, the existence and use of an explicit transformation that 
lowers the cutoff allows one to directly study scaling behavior before using the hamiltonian 
to fit data. 



When the SRG is used with coupling coherence [112, 113], which we explain in the text, 
it allows one to construct effective theories with the same number of free parameters as the 
underlying 'fundamental' theory. For the delta-function examples there is one fundamental 
parameter, the strength of the regulated delta-function as the cutoff is removed. In EFT this 
single parameter is replaced by an arbitrarily large number of effective parameters that are 
fixed by data. In the SRG with coupling coherence, there is still one fundamental coupling 
and all new couplings are fixed perturbative functions of the fundamental coupling. It is 
the renormalization group flow of the added couplings, and a boundary condition that they 
vanish when the fundamental coupling is taken to zero, that fixes their dependence on the 
fundamental coupling. 

The examples we use in this article do not illustrate non-Gaussian fixed points, so their 
scaling properties are driven by naive dimensional analysis. However, we will see that even in 
these cases scaling behavior of effective hamiltonians derived using a perturbative similarity 
renormalization group can be very complicated. We will see that errors in EFT scale as 
inverse powers of the cutoff, which is what one expects since these operators are irrelevant. 
However, in the perturbative SRG there are additional errors arising from the approximate 
treatment of the fundamental running coupling and the approximate treatment of the relation 
between this coupling and the new couplings of irrelevant operators. 

The easiest way to understand errors in the perturbative SRG is by comparison with EFT. 
The same set of operators are found in the SRG hamiltonian as in the EFT hamiltonian, but 
the coefficients are approximated using the perturbative SRG equations instead of being fit 
to data. 

In a realistic calculation the marginal coupling, which corresponds to the strength of the 
regulated delta function, would be fit to data. In order to clearly illustrate the logarithmic 
errors that result from using the perturbative SRG equations, we approximate this marginal 
coupling in this article rather than renormalizing it nonperturbatively by fitting data. The 
strengths of the irrelevant operators, which correspond to derivatives of the regulated delta 
function, are approximated using expansions in powers of the approximate running coupling 
that are fixed by coupling coherence. The approximate running coupling differs from the 
exact running coupling by inverse powers of logarithms of the cutoff, and the error analysis 
for the binding energy displays the resultant inverse logarithmic errors in addition to power- 
law errors similar to those seen in EFT. In addition there are errors in the strengths of the 
irrelevant operators resulting from using a truncated expansion in powers of the running 
coupling and an approximate running coupling, both of which introduce inverse logarithmic 
errors in addition to the power-law errors seen for these operators in EFT. 

We do not argue that EFT or that SRG is to be preferred. Both allow us to produce 
effective hamiltonians whose errors can be understood and thereby controlled. In some 



problems it will be clear that it is much easier to produce an approximate hamiltonian 
using EFT, but this requires the calculation of scattering observables and in field theories 
it may require the use of much of the existing data for a system just to fix the parameters. 
In the SRG the number of parameters equals the number of parameters in the underlying 
'fundamental' theory, but the SRG equations become extremely complicated as one goes 
beyond second order. In addition EFT can be used even when there is no underlying theory 
with the same degrees of freedom, but coupling coherence should fail in this case. This issue 
deserves investigation. 

There are interesting questions about whether the SRG formalism can be used with 
coupling coherence when there is no clear connection between the effective theory and any 
underlying fundamental theory. For example, the number of couplings in chiral perturbation 
theory far exceeds the number of parameters in QCD. If chiral effective theories are exact 
representations of QCD in the low energy limit, coupling coherence would lead one to expect 
that the number of independent couplings may be the same, so that there must be constraints 
on the coupling in chiral field theories that have not been determined. This article does not 
explore some of the most interesting questions, but it is intended to clearly lay the groundwork 
required to initiate such investigations. 

2 Effective Field Theory Approach 

EFT provides a systematic procedure for improving the approximate description of composite 
effective degrees of freedom and their dynamics by increasing the number of local effective in- 
teractions. The effects of high-energy (short-distance) degrees of freedom are incorporated in 
effective interactions whose strength is adjusted to fit appropriate low-energy (long-distance) 
results. The keys to maintaining predictive power in EFT are approximate locality and an 
expansion in powers of a small ratio of low-to-high energy scales or equivalently, short-to-long 
distance scales. 

We begin by discussing the formulation of an effective description for a system of two 
non-relativistic particles that interact via a Dirac delta-function potential in D-dimensions, 



following Lepage's presentation of effective field theory [37]. As mentioned in the Introduc- 



tion, there is a large literature on delta-function interactions where all of the basic results we 



need can be found [14]-[ 31 1 . The hamiltonian describing the relative motion of the system 
is given by 

H = -V 2 r -a 5^(r). (1) 

We use units in which h = 1 and choose the coupling oq to be dimensionless. In the 
D = 1 case the kinetic energy operator is implicitly divided by 2m = 1, where m is the 
two-particle system reduced mass. In the D = 2 case a mass scale is not necessary, since 
with a dimensionless coupling the hamiltonian is scale invariant. 

First, we define a regularized delta- function potential using a cutoff A that eliminates 
contributions from processes involving momentum transfer |q| > A (short-distance dynam- 
ics): 



The cutoff also regulates the interaction at r = 0, avoiding possible divergences in the 



calculation of matrix elements. By Fourier transforming we obtain 
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where q = p — p' is the transferred momentum. 

Next, we introduce an effective potential that consists of a series of approximately local 
contact interactions 
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In momentum space, the potential consists of the Fourier transform of the regularized 
delta-function and an infinite number of approximately local effective interactions (non- 
renormalizable) 
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The effective interactions correspond to the derivatives of the delta-function in momentum 
space and model the effects of the excluded high-momentum transfer contributions, system- 
atically eliminating the dependence of the observables on the cutoff. 

To simplify the calculations, we drop terms such as p.p', obtaining an equally satisfactory 
effective potential that is easier to work with because it is completely separable: 
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The parameters Cj in the expansion can be determined order by order by fitting data for 
low-energy processes. Here we use as "data" the values for the inverse "on-shell" K-matrix, 
calculated with an exact theory. We follow the method developed by Steele and Furnstahl 
|52j| , determining the parameters by requiring the EFT calculation to give the same result 
as the calculation using the exact theory to a given order of accuracy. Instead of using the 
T-matrix, we iterate the effective potential via the Lippmann-Schwinger equation for the K- 
matrix, which is real and therefore more convenient for numerical calculations of scattering 
observables. 

The Lippmann-Schwinger equation for the K-matrix with the effective potential, gener- 
alized to D dimensions is given by 



K x ( P ,p';k)=V x (p,p')+vJ d D q 



Vx(p,q) 
k 2 — a 2 



K x (q,p';k) 



(7) 



where V denotes the principal value. In the D = 2 case only S-wave scattering occurs, 
so the above equation refers to the S-wave potential (angular variable integrated out) and 
d 2 q = q dq. To obtain the K-matrix we solve the Lippmann-Schwinger equation (non- 
perturbatively), using the technique described in Refs. [31, 47, f)3, 60], which is suitable for 



separable potentials. The effective potential given in Eq. (||) (truncated at order p 4 ) can be 
written in the form 
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where A,,- are the matrix elements of 

/ C (A) C 2 (A)/2A 2 C 4 (A)/4A 4 \ 

C 2 (A)/2A 2 C^(A)/2A 4 

\ C 4 (A)/4A 4 / 

The solution of the Lippmann-Schwinger equation can then be written in the form 
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The unknown matrix r(k) satisfies the equation 
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Solving Eq. (11) for r analytically and substituting the solution in Eq. ( |i0|) we obtain the 
effective K-matrix. 

Including only the first operator (with Co) in the effective potential we obtain 



#i°W;fc) 



C (A) 



e 2\ 2 e 2\ 2 



1 - Co (A) I (k) 
Including the first two operators (with Cq and C2) we obtain 
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We then fit the difference between the effective and exact inverse "on-shell" K-matrix 
(p 2 = p' 2 = k 2 ) to an interpolating polynomial in k 2 /X 2 to highest possible order, 
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(16) 



The coefficients A^ are minimized with respect to the variation in the parameters Cj of the 
effective potential. The number of coefficients that can be minimized is given by the number 
of parameters appearing in the effective potential. This method is robust and numerically 
more stable than the matching at discrete momenta used by Lepage, allowing extension to 



higher orders [52]. 



By adjusting only the parameter Co (A) we should eliminate the leading error. As each 
term is added to the effective potential, followed by the adjustment of the respective pa- 
rameter, we expect the errors in the "on-shell" K-matrix to be systematically reduced by 
powers of k 2 /X 2 . Thus, in a log-log plot for A[1/.RT] we expect to obtain straight lines with 
slope given by the dominant power of k/X in the error. We know of no rigorous proof of this 
expectation in EFT, although it may be possible to show that each Ai can be taken to zero 
with a complete set of operators to the appropriate order in p 2 jX 2 in the hamiltonian. Note, 
however, that the number of couplings, Cj, grows more rapidly than the number of A[s and 
the relation between them is nonlinear. We assume that we can adjust operators of order 
(p 2 /X 2 ) m to remove errors in the inverse "on-shell" K-matrix of order (k 2 /X 2 ) m , and all of 
our results confirm this. 

The expansion given by Eq.( |l6]), at best, becomes invalid when the momenta involved 
are of the same order as the cutoff, where one expects the short distance effects to be directly 
resolved. This point corresponds to the radius of convergence of the effective theory. 

The identification of this power counting for the order of errors implies that the same sys- 
tematic dependence on momentum is expected in the evaluation of errors for other low-energy 
observables. It is important to observe that this power counting applies to the predictions 
for the observables and not to the potential. As new terms are added and the corresponding 
couplings are adjusted, the lower-order couplings change. This occurs because the trun- 
cated potential still contains all orders in (p/X) 2 due to the regulating function, so when 
the Lippmann-Schwinger equation is solved non-perturbatively all powers of momentum are 
generated and there is no simple linear relationship between powers of momentum in the 
potential and powers in the scattering observables. 

One of our goals is to compare the EFT and the SRG methods. In principle, this could be 
achieved by comparing the results of the error analysis for low-energy scattering observables 
obtained in each case. However, the calculation of the K-matrix (or T-matrix) with the 
renormalized potential obtained using Wegner's transformation |8(| requires the numerical 
solution of an integral equation involving a non-separable kernel, which is unstable when 
simple techniques are applied. Since we expect the bound-state energy and the low-energy 
scattering observables to have similar error scaling (after fixing the effective potential to a 
given order), we focus on the bound state problem, which is numerically much more stable 
and allows a satisfactory error analysis. Using the effective potential (truncated at a given 
order) with the parameters fixed by fitting the inverse "on-shell" K-matrix, we diagonalize the 
effective hamiltonian numerically, obtaining the bound-state energies. This is implemented 
by discretizing the Schrodinger equation. 

3 Similarity Renormalization Group 

In this section we review the general formulation of the SRG developed by Glazek and Wilson 



1 78, 7£] and a specific transformation developed by Wegner pi. The reader may wish to 



skip the general formulation on a first reading. 

3.1 Glazek- Wilson Formulation 

Consider a system described by a hamiltonian written in the form 

H = h + V, (17) 

where h is the free hamiltonian and V is an interaction. 

In general, the hamiltonian can couple states of all energy scales and such couplings can be 
a source of ultraviolet divergences. The goal of the SRG is to obtain an effective hamiltonian 
in which the couplings between high and low-energy states are removed, while avoiding 
any problems from small energy denominators in effective interactions. The procedure is 
implemented by a unitary transformation that generates effective interactions that reproduce 
the effects of the eliminated couplings. The effective hamiltonian cannot produce ultraviolet 
divergences at any order in perturbation theory as long as its matrix elements are finite. 

In our discussion we will use the basis of eigenstates of the free hamiltonian, 

h\i >= €i\i > . (18) 

We start by defining a bare hamiltonian, H\, regulated by a very large cutoff A (here with 
dimensions of energy) on the change in free energy at the interaction vertices, 

H A = h + V A , (19) 

Va = /aF a , (20) 

F A = v + Hf, (21) 

where f A is a "similarity function", V A is defined as the reduced interaction and H^ are 
counterterms that must be determined through the process of renormalization in order to 
remove A dependence in physical quantities. 

The similarity function f A regulates the hamiltonian by suppressing matrix elements 
between free states with significantly large energy difference and acts in the following way: 

< ilfAlJ > = e i s a + /a(^ - e i) < i\V A \j > 

= tiSij + fMiVw (22) 

Typically, the similarity function is chosen to be a smooth function satisfying 

(0/a( £ j - e j) -*■ 1) wnen \ e i ~ e j\ « A , 

(ii)f A (ei -ej)-* 0, when |e» - e,-| >> A . (23) 

In several papers the similarity function has been chosen to be a step function. Although 
useful for doing analytic calculations, such a choice can lead to pathologies [87]. 



The similarity transformation is defined to act on the bare regulated Hamiltonian, Ha_, 
lowering the cutoff down to a scale A: 

H x = U(X,A)H A UH\,A). (24) 

10 



The renormalized Hamiltonian can be written in the general form 

H x = h + V x , (25) 

V x = fxV x . (26) 

The transformation is unitary, so H\ and H\ produce the same spectra for observables. Also, 
if an exact transformation is implemented, the physical predictions using the renormalized 
Hamiltonian must be independent of the cutoff A and Hff is chosen so that they also become 
independent of A as A — ► oo. 

The unitarity condition is given by: 

U(X,A)UHX,A) = U\X,A)U(X,A) = l. (27) 

The similarity transformation U can be defined in terms of an anti-hermitian operator T\ 
(TJ = —T\) which generates infinitesimal changes of the cutoff energy scale, 

U(X,A) = TexpU T x > dX'\ , (28) 

where T orders operators from left to right in order of increasing energy scale A'. Using 

rA = y(A , A) ^> = _^V (A , A) , (29) 



and the unitarity condition Eq. (27), we can write Eq. (24) in a differential form, 

dH x 



dX 



[Hx,T x ] . (30) 



This is a first-order differential equation, which is solved with the boundary condition 
H\\ x _ tA = H\. The bare Hamiltonian is typically given by the canonical Hamiltonian plus 
counterterms that must be uniquely fixed to complete the renormalization. 

The operator T\ is defined by specifying how Vx and h depend on the cutoff scale A. For 
simplicity in this paper, we demand that h is independent of A, although this may not lead 
to an increasingly diagonal effective hamiltonian in all cases. We also demand that no small 
energy denominators can appear in the hamiltonian. These constraints are implemented by 
the conditions 

dh . , 

Tx s °> (31) 

^ = [Vx,T x ]. (32) 

To obtain the renormalized Hamiltonian perturbatively, expand 

(33) 
(34) 
(35) 



_ _(1) _(2) 

Vx = V x +V X +•• 


5 
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H c ; = H ^ ct + H ^ ct 
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where the superscripts denote the order in the original interaction, V. A general form of 
these effective interactions is 

Vx = E <W) / dX' WJP^+H?-* , (36) 

j,k=i Jx 

— f 1 ) 

for z = 2, 3, • • -, with V^ = v. For instance, the explicit form of the second-order effective 

_(2) 

interaction V x is 

/ (AA) (AA)\ 



ifc ~jfc 



where 



(AA) _ /" A „/ , d/A'fci , Q5A 

ft* = ^ ^ / Vifc _ , (38) 

Ay = ej - ej . (39) 

The counterterms i?A can be determined order-by-order using the idea of coupling- 



coherence [112, 113]. This is implemented by requiring the hamiltonian to reproduce itself 
in form under the similarity transformation, the only change being explicit dependence on 
the running cutoff in the operators and the implicit cutoff dependence in a finite number 
of independent running couplings. All other couplings depend on the cutoff only through 
their dependence on the independent couplings. In general, we also demand the dependent 
couplings to vanish when the independent couplings are taken to zero; i.e, the interactions are 
turned off. If the only independent coupling in the theory is a\, the renormalized hamiltonian 
can be written as an expansion in powers of this coupling: 

H x = h + a A oW + a\0^ + ... . (40) 

In this way, the effective hamiltonian obtained using the similarity transformation is com- 
pletely determined by the underlying theory. The procedure can be extended to arbitrarily 
high orders, although it becomes increasingly complex both analytically and numerically. 

3.2 Wegner Formulation 



The Wegner formulation of the SRG [ BO [ is defined in a very elegant way in terms of a flow 

he Glazek- Wilson formalism 173, 173] , 

[H S ,T S \. (41) 



equation analogous to the SRG Equation in the Glazek- Wilson formalism 173, Yt 



ds 

Here the hamiltonian H s = h + v s evolves with a flow parameter s that ranges from to oo. 
The flow-parameter has dimensions 1/ (energy) 2 and is given in terms of the similarity cutoff 
A by s = 1/A 2 . 

In Wegner's scheme the similarity transformation in defined by an explicit form for the 
generator of the similarity transformation, T s = [H s ,Hq], which corresponds to the choice 
of a gaussian similarity function with uniform width. In the original formulation, Wegner 
advocates the inclusion of the full diagonal part of the hamiltonian at scale s in Hq. For a 
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perturbative calculation of H s , we can use the free hamiltonian, Hq = h. With this choice, 
the flow equation for the hamiltonian is given by 

dH 9 



ds 



[H 8 ,[H s ,h]]. (42) 



The reduced interaction, V s ij (the interaction with the gaussian similarity function fac- 
tored out) is defined by 

V sij = f sij V sij , (43) 

f sij = e~ s ^ • (44) 

Assuming that the free hamiltonian is independent of s, we obtain the flow equation for the 
reduced interaction, 

^ = £ (A lk + A jk ) V slk V skj e ~ 2sA ^ , (45) 

as k 

where we use A?- — A 2 k — A 2 - k = —2A ik Aj k . We should emphasize that this is an exact 
equation. 

To solve this equation we impose a boundary condition, H s \ s ^ g = H So . Then, we make 

a perturbative expansion, 

— — « — ( 2 > , x 

V S = V S +V S +■■■ , (46) 

where the superscript implies the order in the bare interaction V s . It is important to 

observe that counterterms are implicit in the bare interaction and can be determined in the 

renormalization process using coupling coherence. 

At first order we have 

_ (i) 
dV... 

—T 2 - = , (47) 

as 

which implies 

V% = V SQij , (48) 

where s is the final scale. Because of the dimensions of the flow parameter we have s > s , 
corresponding to a smaller cutoff. The "no cutoff limit" corresponds to s — ► 0. 
At second order we have 

_(2) 

—^ - Y, (\k + A jk ) V SQtk V SQkj e- 2sA ^« . (49) 



ds k 



Integrating, we obtain 



—(2) 1 V^— — ( l l \ 



2s A xk A jk _ -2sA ik A jk 



-"^u—* ik —j re p 



(50) 



By construction, the Wegner transformation is unitary and avoids small energy denom- 
inators. The Wegner transformation is one of the Glazek- Wilson transformations, with the 

—A 2 /A 2 
similarity function chosen to be /... = e «' 
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3.3 Strategy 

In our applications of the SRG we use Wegner's transformation. The renormalized hamilto- 
nian for the non-relativistic delta-function potential in D-dimensions is given by 

H x (p,p')=p 2 5 ( - D \p-p') + e~ iE ^ 1]1 [F A (1) (p,p') + y A (2) (p, P , ) + -] , (51) 

where 






^W) = -7^. (52) 



vfW) = al t FP(p,p'), (53) 

vfW') = al t F^(p,p'). (54) 

Here A is a momentum cutoff (as opposed to the energy cutoff discussed above) related to 
the flow parameter by s = 1/A 4 . The index i denotes the order of the calculation for the 
running coupling. 

The renormalized hamiltonian can be used to compute eigenvalues and eigenstates. Since 
the hamiltonian is derived perturbatively we expect cutoff dependent errors in the observ- 
ables. Formally, we can regroup the terms in the renormalized hamiltonian and write it 
as a momentum expansion similar to the one for EFT in Eq. (J6|). The difference is that 
the expansion parameters are analytic functions of the running coupling a\. Expanding the 
operators Fg (p, p') in powers ofp 2 /X 2 we obtain 

F<»>(P,P<) = .„ + , 2 <^> + , <^tp + « !g£ + . . . , (55) 



where the Zi s are constants. Regrouping the terms we obtain 

Hx(p,p) = P 2 <5 (D) ( P - p') + e~^^ 
where 



(p 2 +p' 2 ) 
9o{a>\) + #2(«a) ^ h • 



(56) 



gi(ot\) = Oi(A) a x + 6, (A) a\ + ... . (57) 

For D = 2 all aj(A), bi(A), etc. are simply constants, while for D = 1 they are functions of 
A because of the implicit mass scale which is hidden when we set 2m = 1. 

We can identify three interdependent sources of errors in the perturbative similarity 
renormalization group when the hamiltonian given by Eq. (|56| ) is truncated and used to 
compute a physical quantity: 

a) errors introduced by the truncation of the hamiltonian at a given order in p 2 /X 2 , which 
are analogous to the errors in EFT; 

b) errors introduced by the truncation of the hamiltonian at a given order in the running 
coupling a\ ; i, which correspond to the use of an approximation for the functions gf, 

c) errors introduced by the approximation for the running coupling ot\^. 
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In the actual calculation using the hamiltonian given by Eq. ( |5l| ) errors of type (a) do not 
appear directly because we do not truncate the operators that appear in the hamiltonian. 
However, errors of type (b) can be understood as coming from approximating the couplings in 
front of the operators in Eq. (|56|). In EFT these couplings are fit to data, but we approximate 
them as in Eq. ( p7|) and this leads to an approximate cancellation of the power-law errors 
completely removed by these operators in EFT. Errors of type (c) appear in our calculations 
only because we do not fit the canonical coupling to data at each scale, but fix it at a given 
scale and evolve it perturbatively from that scale. The strategy we would use for a realistic 
theory (e.g., QED and QCD) is the following: 

1) Obtain the renormalized hamiltonian using the similarity transformation and coupling- 
coherence, truncating the hamiltonian at a given order in powers of a\^. 

2) Fix the coupling a\ by fitting an observable (e.g., a bound-state energy). 

3) Evaluate other observables (e.g., scattering phase shifts). 

As pointed out before, the evaluation of scattering observables with the similarity hamil- 
tonian with standard techniques is complicated and so in our examples we focus on the bound 
state errors. We fix the coupling at some scale using a given renormalization prescription and 
use the flow-equation to obtain the coupling as a function of the cutoff A to a given order. 
We then perform a sequence of bound-state calculations with better approximations for the 
hamiltonian such that the errors in the bound-state energy are systematically reduced. Once 
the sources of errors are identified, it becomes relatively simple to analyze order-by-order 
how such errors scale with the cutoff A. In principle, to completely eliminate the errors 
proportional to some power m in the momentum expansion we should use the similarity 
hamiltonian with the exact running coupling (renormalized to all orders) and include the 
contributions up to 0(p m /X m ) coming from all effective interactions (all orders in a\). Since 
we use a perturbative approximation for the running coupling and only a finite number of 
effective interactions, we do not expect the scaling analysis for the errors in SRG to be as 
simple as in the EFT approach. Some details of this scaling analysis are presented later for 
the specific examples we work out. We should emphasize that in a realistic calculation we 
would fit the coupling a\ to an observable. This nonperturbative renormalization eliminates 
the dominant source of errors we display in SRG calculations in this paper. We choose to 
renormalize the coupling perturbatively in this paper because the only observable we com- 
pute is the single bound state energy of a delta-function potential, and fitting this energy 
would prevent us from displaying errors. 

4 One-Dimensional Delta-function Potential 

Consider a system of two non-relativistic particles in one dimension interacting via an at- 
tractive Dirac delta-function, contact or zero range potential. The Schrodinger equation in 
position space (with Ira = 1 and fi = 1), separating out the center of mass motion, is given 
by 

/L',T,," .'I 

a 5 (1) (x) *(x) =EV(x) . (58) 



d 2 ^(x) 



dx 2 
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The coupling constant oq completely fixes the underlying theory. 

4.1 Exact Solution 

The exact solutions of this equation are familiar to all students of quantum mechanics. 

a 2 

There is one bound state, (E < 0) with binding energy Eq = —E = -£■ and normalized 
wave-function 

*(x) = J^-e- S £W. (59) 

In the scattering problem (E = k 2 > 0) the Schrodinger equation is equivalent to a Lippmann- 
Schwinger equation, whose solution with outgoing waves is given by 



* k (x)= '"' 



2tt 



2ik 



2ik + «o 



(60) 



The problem is easily solved in momentum space, also. We do this to lay the ground- 
work for our first illustration of renormalization techniques. The corresponding Schrodinger 
equation in momentum space is given by 

p 2 $(p) - £° f°° dq <%) = E $(p) , (61) 

27T J-oo 

where $>(p) is the Fourier transform of the position space wave- function, 

$(p) = -£= r dx 9{x) e~ ipx . (62) 

V 27T J-oo 

Identifying the integral in Eq. ( |6l|) as the position space wave-function at the origin 

1 Z" 00 
-= / dq <%) = *(0) , (63) 

V ZTT J — oo 

we obtain 

p 2 $(p) - -g= tt(0) = E $(p) . (64) 

V27T 



First, we consider the bound-state problem. By rearranging the terms, we obtain 

_ _a^_ tt(Q) 
* {P) - V2^(P 2 + E )- (65) 



Integrating the above equation in p, we obtain 



2vr i-oo P (p 2 + E ) 



1 = S= / ^7-^^T- (66) 



2 

This integral fixes the binding energy, Eq = -£-. The corresponding bound-state wave- 



function in momentum space is 



1 9 FA 

W v^(P 2 + ^o) { ' 
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Solving for the scattering states we obtain the solution corresponding to outgoing waves 



<Mp) = <s(p-*0 + g 



2ik 



2ik + an 



1 



(p 2 — k 2 — ie) 



(68) 



where k = yE. 

All scattering observables can be obtained from the T-matrix, which is given by a 
Lippmann-Schwinger equation: 



/oo 
dq 
-oo 



V(p,q) 



k 2 



+ ie 



T(q,p';k) 



The "on-shell" (p = p = k ) T-matrix for a Dirac delta- function potential is 



Tn(k) 



2n 



2ik 



2ik + an 



(69) 



(70) 



A more convenient way to treat the scattering problem for our purposes is to use the 
K-matrix, which satisfies the Lippmann-Schwinger equation 



K(p,p';k) =V(p,p')+V 



dq YML K ^ p '. k ) 

k z — o z 



(71) 



This equation is similar to that for the T-matrix, except standing-wave boundary conditions 
are imposed so that the usual ie prescription is replaced by the principal value V . As a 
result, the K-matrix is real. The on-shell K-matrix is given by 



The "on-shell" K-matrix and T-matrix are related by the expression 

T o(k) 



(72) 

(73) 

(74) 
(75) 

(76) 



Although the one-dimensional problem can be solved exactly without regularization and 
renormalization, for the sake of further comparison with the two-dimensional case and to il- 
lustrate the SRG and EFT techniques, it is worth analyzing what happens when a momentum 
cutoff is introduced. If the integral in Eq. (pq) is regulated by a cutoff A, 



(77) 









J1 uw — 1 


fn(k) ■ 


Using either 






k cot 6o(k) — ik 


k 2 i 

vr T (k) 


or 






k cot 8n(k) = - 


k 2 1 




vr K (k) ' 


the phase-shifts 


are 


easily 


calculated: 

cot 8n(k) 


2k 
an ' 



2vr La P (p 2 + E ) ' 
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we obtain a transcendental equation, whose solution gives a cutoff-dependent binding energy: 

A \ 



7T 



Eq(A) = arctg 



(78) 



Clearly, the cutoff dependence of the bound state energy is eliminated if we let the coupling 
run like 



"0 — ► «A 



arct e v ./Eb 



(79) 



For A >> -\/Eq this cutoff dependence can be eliminated perturbatively by using the expan- 
sion of «a in powers of ^/Eq/A to a given order: 



a A = 2VE 



1 + l^ + ±Sl + o(^ 



7T A 7T 2 A 2 

In terms of the "bare" coupling, qq, we obtain 



(80) 



«A = "0 



1 Qo 1 Oq 



n; 



7T A 7T 2 A 2 \ A 3 



(81) 



Of course a a — ¥ &0 when A — > oo. 



4.2 Effective Field Theory Method 

We now apply the EFT approach to the one-dimensional delta-function problem. The power 
counting scheme in this case is very simple and can be understood easily by analyzing how the 
errors scale with the cutoff A. For simplicity, we consider the errors in the inverse K-matrix. 
As described before, the EFT method is designed to systematically eliminate errors order 
by order in k/X. Beginning with the leading-order prediction, we use the effective potential 



with one parameter. Evaluating the difference A [l/li^ '] = 1/K^ 
it in powers of k/X we obtain: 



(()) 



1/Kq and expanding 



A 



K(°) 



1 2^F | 2tt\ / 1 
Co A a / V C 



TV 



3A 



k 2 „{k 4 " 
X 2+ °W 



(82) 



Note that the regulator generates only power-law errors, which indicates that there are no 
divergences as A — ► oo and renormalization is required only to remove finite errors. Clearly, 
the leading order error is eliminated by choosing 



C n 



— ao/2ir 
1 ,R A 



(83) 



and the remaining errors are dominated by the term 



O 



A 2 



2tt 4^ 
uq 3A 



A 2 • 



(84) 
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For the next-to-leading order prediction, we use the effective potential with two parame- 
ters, obtaining the expansion 



A 



1 



KM 



k 2 I A- 4 \ 

Ao + Atjp + ol-g , (85) 



where 



A = - 2 ( l + ^) +2 ^7^ ^^ + - (80) 

-2 Co + %^ «o 



8C 2 ^ + 4C A(3 + 3^ + ^ 



3 A -2C + 



4 A 



r v / 1 9 i 33 C 2 y^F , 15 C| 7T 5 Cf tt 3 / 2 
G 2 A ^ 12 H 2A h 2A^ H 4A3— 

3A (-2Co + %^) 2 



(87) 



The parameters Cq and C2 are determined by solving the system of two coupled non-linear 
equations given by 

A) = 0, A 2 = 0. (88) 

One of the solutions gives a parameter Cq that diverges when A — ► 00 and so is not acceptable. 
The other solution, 



Gi 



-270F a A 2 + 16 a 2 Q A + 12 A 3 tt 
2^F(4 a 2 , - 9 a A 0r + 6 A 2 tt) 



^ (ao-A y^F) 2 (4 a 2 _ 9 ao A 0F + 6 A 2 0F) 

4 a 2 - 9 a A 0r + 6 A 2 vr ' ^ ' 

-8 a§ v^ A + 18 a vr A 2 - 12 A 3 tt 3 / 2 
7r(4 a^ - 9 a A ^/7r + 6 A 2 7r) 



C 2 



2^/24 a 4 A 2 tt - 102 a 3 A 3 vr 3 / 2 + 168 a 2 A 2 tt 2 - 126 a A 5 vr 5 / 2 + 36 A 6 tt 3 
+ vr(4 ag - 9 a A V^ + 6 A 2 vr) ' 

(90) 

is such that Cq — ► — ao/2ir, C2 — ► when A — ► 00 and leads to the elimination of errors up 
toO(g). 

This procedure can be systematically extended to higher orders, although with expo- 
nentially increasing algebraic complexity. It is important to observe that the scaling of the 
parameters, and hence the power-counting scheme, depends on the regularization procedure. 
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Figure 2: The leading-order and next-to- leading-order EFT expansion parameters for the one- 
dimensional delta-function potential. 



In Fig. 2 we plot the parameters obtained in the leading order and next-to-leading order 
calculation as a function of the cutoff. Note that when we add the term proportional to Ci 
the parameter Cq changes. The numerical evaluation here and in what follows is implemented 
with the underlying D = 1 theory fixed by choosing ao = 2, which corresponds to an exact 
binding energy Eq = 1. 

We proceed by evaluating the errors in the binding energy, using the effective potential 
with one, two and four parameters adjusted to fit the inverse K-matrix. Instead of solving 
systems of non-linear coupled equations, we fix the coefficients by applying the minimization 
procedure described in Section 2, which is simpler numerically and leads to the same results. 

We solve the Schrodinger equation numerically for different values of the cutoff A. First, 
we introduce an integration cutoff Aq >> A, such that 



P $\(P) + / dq V x (p,p') $ A (ff) = E x <S>\(p) . 

J -An 



We discretize the Schrodinger equation (choosing a gaussian quadrature) 



N 
3=1 



PiPj Sij \-V\(pi,pj) 



IV 



$x(Pj) = Ex $x(Pi) , 



obtaining a matrix equation 



(M - E x 1)* A = 



(91) 

(92) 
(93) 
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Figure 3: The EFT error in the binding energy for the one- dimensional delta-function using one, 
two and four parameters. The exact theory is fixed by choosing a = 2. 



where 



M, 



PiPj <% — + V x (pi,pj] 



u< 



(94) 



Next, the matrix M is symmetrized, reduced to tridiagonal form and the resulting matrix 
equation is solved by the QL method. The smallest (negative) eigenvalue corresponds to 



the bound-state energy. We then refine the result using Richardson extrapolation [114] in 
order to eliminate the dominant errors due to the integration cutoff, which are proportional 
to Eq/Aq. Removing these errors should not be confused with removing A-dependent errors. 
Finally, we compare the resulting bound-state energy, E\, with the exact value. 

In Fig. 3 we show in a log-log plot the absolute values for the relative errors in the 
bound-state energy as functions of Eq/X 2 . As expected, we obtain straight lines with slope 
given by the dominant power of Eq/X 2 in the error, in agreement with the scaling analysis. 

This result clearly demonstrates the systematic improvement in the accuracy of the ef- 
fective potential for describing the exact theory in the low-energy regime. We also note that 
the errors corresponding to each approximation become comparable as A 2 approaches the 
binding energy, indicating the breakdown of the effective theory. In this way, the error plot 
provides a graphical illustration of the radius of convergence, the point where all the error 
lines converge. 
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4.3 Similarity Renormalization Group Approach 

In the one-dimensional case the canonical hamiltonian in momentum space with a delta- 
function potential can be written as H(p,p') = h(p,p') + V(p,p'), where h(p,p') = p 2 5^(p — 
p') corresponds to the free hamiltonian and V(p,p') = —ao/2ir corresponds to the Fourier 
transform of the delta-function potential. This hamiltonian does not produce divergences so 
an initial cutoff is not required to define the problem. 

The flow equation obtained with Wegner's transformation in terms of matrix elements in 
the basis of free states is given by 



dV s (p,p') 
ds 



(p 2 -p ,2 ) 2 V s (p,p') 



dk (2k 2 -p 2 -p' 2 ) V s {p,k) V s {k,p') , (95) 



with the initial condition H s= o(p,p') = H(p,p'). The reduced interaction V s (p,p') is defined 
such that 

V s (p,p') = e- s ^- p ^ 2 V s (p,p'). (96) 

Now, assuming that h is cutoff independent we obtain the flow equation for the reduced 
interaction, 



dVs_ 
ds 



-2s p 2 p L 



dk{ 2k 2 -p 2 -p' 2 )e- 2s ^- k{ P 2+ P'- 
xV s ( P ,k)V s (k,p'). 
In order to solve this equation, we employ a perturbative expansion 

v.s(p,p) = y, (1 W) + v s ^(p,p') + ■■■ , 

starting with 



v, (1 W) 



"2vr 



(97) 



(98) 



(99) 



Using coupling-coherence [113] we assume a solution in the form of an expansion in powers 
of a s /2ir, 

OO / s Tl 



f,M= ^ + SSJ F * n)M) 



(100) 



with the additional condition that the operators Fs (p,p') vanish when p = p' = or, 
equivalently, when s = 0. Using the solution Eq. ( 100P in Eq. (^) we obtain: 



dVs_ 
ds 



1 da s ^ 1 

2^ ~ds~ + 2-L (2vr)" 

n=2 v ' 



n a 



"" s F^{p,p')+a n s dFs } P ' P ' ] 



da s 
~ds~ ls 



ds 



" dk (2k 2 -p 2 -p' 2 ) e ~MP 2 P' 2 +^-kH P H P ' 2 )] 



OO 

a.s 

~2lT 



+ E 



x 'a^ n 



n=2 



2?r 



F^{p,k) 



a s 
2ir 



+ E 



00 'a, sm 



m=2 



2n 



F^ m \k,p') 



(101) 



This equation is solved iteratively order by order in a s /2n. If a s /2tt is small the operator 
V} (jPiP 1 ) can be identified as the dominant term in the expansion of V s (p,p') in powers of 
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p and p'. In fact, since in the limit when p and p' approach zero all Fg (p,p') vanish, the 
only term left is — a s /(27r), which in the D = 1 case corresponds to a relevant operator (even 
though the coupling is dimensionless, it is multiplied by 2m = 1). The higher-order terms 
correspond to irrelevant operators. Thus, at each order n the solution can be viewed as a 
double expansion obtained by solving the equations for a s and Fg (p,p') separately. 
At second-order we have 



1 da s 



2vr ds (2vr) 



1 2 dFP(p,p>) 



2 "s 



ds 



a. 



li 2) (p,p) , 



where 



4 2 W) 



(2*)* 



dk {2k 2 -p 2 -p' 2 ) e -2 S [pV 2 +fe 4 -fe 2 ( P 2 +^)] 



The equation for the running coupling, a s , is obtained by taking the limit (p,p') — > 0: 



J_^£ = a 2 j( 2) (00) 

2ix ds s s v ' ' 



where 



4 2) (o,o) 



(2*)= 



ci£;2/c 



2 -2sfc 4 



i r(3/4) i 



(2vr) 2 2 3 /4 s 3/4 



Integrating Eq. (|104|) we obtain: 



« s 



«s,2 



'/2 



" i - m ^ (<v< - s y 4 ) ' 

4T(3/4) 

2 3 / 4 ' 



(102) 
(103) 

(104) 
(105) 

(106) 
(107) 



The extra subscript on a displays the order of the right-hand side of Eq. (104). Thus, 
knowing the value of a So for a given so, Eq. ( |106| ) can be used to obtain the running 
coupling a s for any value s. The choice so = corresponds to a SQ = ao ("bare coupling"), 
such that 

a s ,2 



«o 



i-m% a V4 



In terms of the cutoff A we obtain: 



(2), 



«A.2 



a 



!-»»£* 



The equation for Fg (p,p') is given by: 



,(2), 



^ZLMI^m-Vm 



(108) 
(109) 

(110) 



Integrating we obtain: 



ifW) = (2vr) 2 / ds' 1^(0,0) -I^(p,p') 
Jo L 



dfc 






+ 



-2s(k 2 -p 2 ){k 2 -p' 2 ) 



1 



+ 



2{k 2 -p 2 ) 2{k 2 -p' 2 ) 



'MD 
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The integral over k must be performed numerically. 
At third-order we have: 



1 da a 1 n2 dFP(p,p') 

2vr ds (2vr) 



OL 
2 s 



ds 



+ 



da 



(M- 



a, s 



^i 2 W) + 7^« 



1 3 dFi 3) (p,p') 



r/.s 



(27T) 



3 s 



-a. 



lP(p,p') + all^(p,p'), 



where 



4 3) (p,p') 



(2tt) 3 



dk {2 k 2 -p 2 -p' 2 ) e -2s{p> P » + kl-kH P * +P ''- 



FP(p,k)+FP(k,p') 



Taking the limit p,p' — > we obtain: 



1 da s 
27r ds 



«; 



if (0,0)- a* if (0,0). 



where 



ds 



(112) 



(113) 



(114) 



If (0,0) 



(2vr) s 



difc 2/fe 2 e- 2sfc4 



Ff (0,&)+if)(A;,0)' 



For dimensional reasons Eq. ( |114 ) takes the form 






#2 2 



s 3/4 



53 a 3 



,1/2 



(2) (31 

where #2 & n d -B3 can be obtained by evaluating Ig (0,0) and I s (0,0) for s 



(115) 



(116) 



1. In terms 



of the cutoff A we obtain: 



da\ _ 4 
~dA~ = A 2 " 



5 2 al 



A 3 



5 3 a!. 



(117) 



Eq. (114) leads to a transcendental equation which is solved numerically in order to obtain 
the running coupling a s3 . 

—A 2 /A 2 
A simple scaling analysis shows that the similarity function f x . . = e *•?' introduces 

errors that are expected to scale as inverse powers of A. The leading-order errors are pro- 
portional to \/Eq/X and are expected to be eliminated by using the running coupling renor- 
malized to second-order. The next-to-leading-order errors are proportional to Eq/X 2 and are 
expected to be eliminated by using the running coupling renormalized to third-order and in- 
cluding the second-order effective interaction, since both introduce corrections proportional 
to Eq/X 2 , which therefore enter at the same level. Thus, to show the systematic improve- 
ment in the perturbative approximation in our study of the one-dimensional delta-function 
we compute the running coupling to third-order and the operators to second-order. As in 
the EFT calculation, we choose the exact binding-energy Eq = 1 («o = 2) and calculate the 
bound-state energy in a given approximation by solving the Schrodinger equation numeri- 
cally. In Fig. 4 we show the running coupling renormalized to second and third-order as 
a function of the cutoff A. In Fig. 5 we show in a log-log plot the absolute value of the 
relative errors in the binding-energy as a function of the ratio E^ 1 ' /X 2 using the following 
approximations for the interaction: 
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(a) relevant operator with "bare" coupling «o, 

(b) relevant operator with running coupling renormalized to second-order ( 0^,2) 

VxM) 



1 2 l'2\2 

a A , 2 _1il_pJ_ 
2tt 



(118) 



(119) 



(c) relevant operator plus second-order irrelevant operator with running coupling renormal- 

(2) 
ized to third-order (0^,3,-^ ), 

V\{p,p') = 



«A,2 . / Otlambda,2 V ,-,(2) , /, 

+ ( - I F\ (P,P) 



, 2 /2%2 
(P -V ) 

P 



2vr V 2vr 
(d) relevant operator with running coupling renormalized to third-order (a\ t s), 



V\(p,p') 



a A ,3 ^-f) 1 
2vr 



(120) 



fl211 



(e) relevant operator plus second-order irrelevant operator with running coupling renormal- 

(2) 
ized to third-order, (a\^,F^ ), 



V\{p,p) 



a \,3 , ( a lambda,3\ ,-,(2) / a 



(122) 



With hamiltonian (a) (unrenormalized) the dominant error is proportional y^Eo/A, which 
is the leading-order error introduced by the regulator. With the renormalized hamiltonians 
(b), (c), (d) and (e), for intermediate values of A, we obtain straight lines with slope given 
by the dominant power of Eq/X 2 in the error. With hamiltonian (b) the slope is one, cor- 
responding to an error proportional to Eq/X 2 . With hamiltonians (c) and (d) there is only 
a small shift in the errors; there is no significant change in the slope. With hamiltonian 
(e) the slope is two, corresponding to an error proportional to (Eq/X 2 ) 2 . Analogous to the 
EFT calculation, for small values of A the lines converge, indicating the breakdown of the 
perturbative expansion. It is important to note that for very large values of the cutoff A 
the scaling behavior might change, since higher-order contributions not included in a given 
approximation for the hamiltonian can become important. In our calculations we did not 
analyze this regime because numerical accuracy limits the maximum range of cutoffs on the 
integrals. 
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Figure 4: The SRG running coupling for the one-dimensional delta-function. 
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Figure 5: The SRG error in the binding energy for the one-dimensional delta-function using various 
approximations for the SRG hamiltonian. The exact theory is fixed by choosing ao = 2. 
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5 Two- Dimensional Delta-function Potential 

We now consider the case of two nonrelativistic particles in two dimensions interacting via 
an attractive Dirac delta-function potential. The Schrodinger equation for relative motion 
in position space (with ft = 1), can be written as: 

- V 2 ^(r) - a 5 (2) (r) *(r) = E tf (r) . (123) 

Both the delta-function potential in two dimensions and the kinetic energy operator scale as 
1/r 2 , therefore, the coupling cto is dimensionless. As a consequence, the hamiltonian is scale 
invariant (i.e., there is no intrinsic energy scale) and we can anticipate the presence of loga- 
rithmic ultraviolet divergences, analogous to those appearing in QED and QCD. The problem 
requires renormalization. In this subsection we present the standard method that produces 



an exact solution analytically, using simple regularization and renormalization schemes [20 1 . 



5.1 Exact Solution 

We start with the Schrodinger equation in momentum space, 

P 2 <&(p) " j^y J d 2 q $(q) = E *(p) , (124) 

where 3>(p) is the Fourier transform of the position space wave- function, 

$(p) = -L / d 2 r tf(r) e~ ipr . (125) 

2-7T J 



As a consequence of scale invariance, if there is any negative energy solution to Eq. (]124| ) 
then it will admit solutions for any E < 0. This corresponds to a continuum of bound 
states with energies extending down to — oo, so the system is not bounded from below. By 
rearranging the terms in the Schrodinger equation we obtain 

, , an ^(0) 

where ^(0) is the position space wave- function at the origin and Eq > is the binding energy. 
To obtain the eigenvalue condition for the binding energy, we can integrate both sides of 
Eq. (|I2§): 

1 = £° f°° dpp . l — . (127) 

2vr Jo (P 2 + E ) V 7 

The integral on the r.h.s. diverges logarithmically, so the problem is ill-defined. 

The conventional way to deal with this problem is renormalization. First, we regulate 
the integral with a momentum cutoff, obtaining 

1 = S C * P 73-f-ETT = £ ln ( 1 + '- I > ( 128 ) 

so that 

Eo = —^ • (129) 

e a c 
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Clearly, if the coupling «o is fixed then Eq — > oo as A — ► oo. In order to eliminate the 
divergence and produce a finite, well-defined bound state we can renormalize the theory by 
demanding that the coupling runs with the cutoff A in such a way that the binding energy 
remains fixed as the cutoff is removed: 



4tt 



ao — ► «a 



In 1 + 



(130) 



The dimensionless renormalized running coupling a\ that characterizes the strength of 
the interaction is therefore replaced by a new (dimensionful) parameter E$ > 0, the binding 
energy of the system. This is a simple example of dimensional transmutation: even though 
the original "bare" hamiltonian is scale invariant, the renormalization procedure leads to a 
scale that characterizes the physical observables. Note that Eq can be chosen arbitrarily, 
fixing the energy scale of the underlying (renormalized) theory. It is also interesting to 
note that the renormalized running coupling a\ vanishes as A — > oo and so the theory is 
asymptotically free. 

This renormalized hamiltonian can be used to compute other observables. The usual 
prescription for the calculations is to obtain the solutions with the cutoff in place and then 
take the limit as the momentum cutoff is removed to oo. If an exact calculation can be im- 
plemented, the final results should be independent of the regularization and renormalization 
schemes. As an example, we calculate the scattering wave function, 



%(p)=^(p-k) + 



«A 



*(0) 



2tt (p 2 — k 2 — i e) 



(131) 



where k = \J E. Integrating both sides over p with a cutoff A in place, we obtain 

-i 



*<°» - i 



47r \ — k A — le 



thus, 

a A *(0) 

In the limit A — > oo we obtain: 



1 

2^ 




4tt 



In 1 + 



A 2 



-i -i 



-k 2 



le 



(132) 



(133) 



«A *(0) 



'-'fi 



l 7T 



The resulting scattering wave function is then given by 



<Mp) = ^ (2) (p-k) + ^- 



2tt (p 2 — k 2 — ie) 



(134) 




n -l 



l IT 



(135) 



It is important to note that only S-wave scattering occurs, corresponding to zero angular 
momentum states. For the higher waves the centrifugal barrier completely screens the delta- 
function potential and the non-zero angular momentum scattering states are free states. 
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The same prescription can be used to evaluate the T-matrix or the K-matrix. For the 
T-matrix, the Lippmann-Schwinger equation with the renormalized potential is given by: 



T(p, p'; k) = V(p, p') + / d\ 2 V{P ^l , T(q, p'; k) . (136) 

J k z — q z + te 



V (p,cQ 
q 2 + it 

Since only S-wave scattering takes place we can integrate out the angular variable, obtaining 



T0-°W; k) = V^ \p,p') +T dqq V J "^ g ? T^ \q,p'- k) , (137) 

Jo K — q + ie 

where 

V^ HP,P') = -^- (138) 

The Lippmann-Schwinger equation for the "on-shell" T-matrix is given by: 

T^\k) =-^~f- T^\k) I" dq q \ (139) 

27r 2tt Jo k — q + le 

Solving this equation and taking the limit A — > oo, we obtain the exact "on-shell" T-matrix: 

T o(k) = - f J , • (140) 

Here and in what follows we drop the superscript and use the subscript to denote the exact 
quantities. 

In the same way, the S-wave Lippmann-Schwinger equation for the K-matrix is given by 

K(p,p';k)=V(p, P ')+vJ o dqq^^K(q,p';k), (141) 

and the exact "on-shell" K-matrix is given by 

*.<*) = -^ ■ (142) 



The "on-shell" K-matrix and T-matrix are related by 

Tp(k) 
1 - f T (k) 



Ko(k) = -4§77T • (143) 



Using either 

kcot5 {k)-ik = ———, (144) 

vr To{k) 

or 

kcot5 (k) = - — -J—, (145) 

7T Ko{k) 



we can obtain the exact phase-shifts: 



1 I k 2 \ 
cot<Jo = -ln — . (146) 

vr \Eo) 
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5.2 Effective Field Theory Approach 



The EFT power counting scheme in the D = 2 case is also simple. Again, we analyze the 
errors in the inverse K-matrix. For the leading-order prediction, we use the effective potential 
with one parameter. Evaluating the difference A [1/iO '] = l/K^ — 1/Kq and expanding 
it in powers of k/X we obtain: 



A 



K(°) 



By choosing 



J__2 + I ln (^_ 

Co 2 2 \Eq 



Co (A) 



+ 



1 
Co 



k 2 „ fk 4 " 
X 2+ °[x 4 



Mi) 



(147) 



(148) 
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the leading logarithmic error is eliminated and the remaining errors are dominated by the 
term 

_ 1 

~ 2 




7 




1 



A 2 ■ 



(149) 



For the next-to-leading order calculation we use the effective potential with two parame- 
ters, obtaining the expansion 



A 



1 



K& 



k 2 fk 4S 



where 



M 



An 



A + 2C 2 + ^j- 1 
2 Co - ^ 2 



In 



En 



4 C n 2 



2 Co (4 + 2 C 2 + ^) + C 2 



3C? 



+ 12C 2 + 3C| + ^) 



(*7 



2 (2 Co - ^ 



(150) 



(151) 



(152) 



The parameters Cq and C 2 are determined by solving the system of two coupled non-linear 
equations given by 

A) = 0, A 2 = 0. (153) 

The acceptable solution is the one that leads to asymptotic freedom, i.e., (Co,C 2 ) — ► as 
A — > oo (for the other solution, Cq — > 8 and C 2 — > in this limit) . The corresponding 
expressions for Cq and C 2 are algebraically very complicated and so we do not present them 
here. In Fig. 6 we plot the parameters obtained in the leading order and next-to-leading 
order calculation as a function of the cutoff. For the numerical evaluation here and in what 
follows we fix the D = 2 theory by choosing the exact binding energy to be En = 1. 

To evaluate the errors in the binding energy we use the effective potential with one, 
two and four parameters adjusted to fit the inverse K-matrix (using the same minimization 
procedure as in the D = 1 case ). 

We solve the Schrodinger equation numerically for different values of the cutoff A. Intro- 
ducing an integration cutoff Aq >> A, 



f A ° 
P $\(P) + dqq V x (p,p') $*(?) = E x $ x (p) . 

Jo 



(154) 
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Figure 6: The leading-order and next-to-leading-order EFT expansion parameters for the two- 
dimensional delta-function potential. 



Discretizing the Schrodinger equation (choosing a gaussian quadrature) 



N 



PiPj 5ij hpj V x (pi,Pj) 



$ x (Pj) = E x 9\(jh) 



we obtain the matrix equation 



where 



Mij 



(M - E x 1)* A = 



Wj PiPj Sij \-pj V\(pi,pj) 



(155) 



(156) 



(157) 



The bound-state energy, E\, is obtained by following the same procedure as for the D = 1 
case and is compared with the exact value. 

In Fig. 7 we show in a log-log plot the absolute values for the relative errors in the bound- 
state energy as a function of Eq/X 2 . As expected, the result is analogous to the D = 1 case: 
straight lines with slope given by the dominant power of Eq/X 2 in the error. Note that the 
dominant logarithmic errors are completely absorbed by adjusting the parameter Co and 
when more parameters are adjusted we obtain power-law improvement. 
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Figure 7: The EFT error in the binding energy for the two-dimensional delta-function using one, 
two and four parameters. The exact theory is fixed by choosing Eq = 1. 



5.3 Similarity Renormalization Group Approach 

In the two-dimensional case the canonical hamiltonian in momentum space with a delta- 
function potential can be written as 



H(p,p>) = h(p,p') + V(p,p>), 



(158) 



where h(p, p') = p 2 5^ 2 ' (p— p') corresponds to the free hamiltonian and V(p, p') = — oo/(2vr) 2 
corresponds to the Fourier transform of the delta-function potential. 

Integrating out the angular variable, the flow equation obtained with Wegner's transfor- 
mation in terms of matrix elements in the basis of free states is given by 



dV s (p,p') 
ds 



V - p' 2 ) 2 V a (p,p') - dkk (2k 2 -p 2 - p' 2 ) V B {p, k) V s (k,p') . (159) 
Jo 



In principle, we can set the boundary condition at s = (no cutoff), i.e, 

H s=0 (p,p') = H(p,p') =p 2 5^(p-p') - ^ . 



(160) 



However, the hamiltonian with no cutoff produces logarithmic divergences, requiring renor- 
malization. As we will see, the boundary condition must be imposed at some other point, 
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leading to dimensional transmutation [115]. The reduced interaction V s (p,p') is denned such 
that 

V s ( P ,p') = e- s ^~ p ^ 2 V s (p, P '). (161) 

Assuming that h is cutoff independent we obtain the flow equation for the reduced interaction, 



ds 



-2s p 2 p'- 



dk k (2k z - p z - p' z ) e 
xV s (p,k)V s (k,p'). 



,/2\ -2s[k 4 -k 2 (p 2 +p' 2 )] 



(162) 



As in the D = 1 case, this equation is solved using a perturbative expansion given by Eq. 
( |98|) , starting with 

VP(p,p') = -^. (163) 



We assume a coupling-coherent solution in the form of an expansion in powers of a s /2n, 
satisfying the constraint that the operators F$ (p,p') vanish when p = p' = 0, 






2tt ^L\2n 

71=2 



(164) 



Note that the expansion parameter is a s /2ir. 

Using the solution Eq. (163) in Eq. ( |162| ) we obtain 



ds 



1 da s ^ 1 



n a 



"-' "" s Fi n W) + < — } p,p>) 



da s 
ds 



-,{n). 



dA; fc (2/t 2 -p 2 -p' 2 ) e -2 S bV 2 +fe*-fe 2 (p 2 + P ' 



S + Efe) ^>(p.*) -S + Efe) *!r*M 



2vr ^-i V27r7 

ra=2 



2vr -H, V2vr 

m=2 



(165) 



This equation is solved iteratively order-by-order in a s /2ir. Again, if a s /2ir is small the 
operator V} (p,p r ) can be identified as the dominant term in the expansion of V s (p,p') in 
powers of p and p' . In the D = 2 case this operator corresponds to a marginal operator 
(since the coupling is dimensionless and there is no implicit mass scale). The higher-order 
terms correspond to irrelevant operators. 
At second-order we have 



J_da 1+ 1 a 2 dF^(p,p') 



2vr ds (2vr) 
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(166) 



where 



4 2 W) 



i 



(2^) 2 Jo 
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dk k (2k 2 -p 2 - P ' 2 ) e -MP 2 P' 2 +^-kH P 2 +P > 2 )] 



2s p 2 p' 2 



4s 



(167) 
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The equation for a s is obtained by taking the limit (p,p') — > 0, 

_1_ dOs_ 2 (2) 

2tt ds ' ' l ' j ' 



where 



If (0,0) 



1 1 



(2vr) 2 4s 



Integrating Eq. (|168|) from sq to s, 



a 



a s ,2 



so 



^ln(^) ' 



In terms of the cutoff A we obtain 



"A,2 



"An 



1 + ^^(A)" 



(168) 
(169) 

(170) 
(171) 



In principle, knowing the value of a so for a given so we can determine the running coupling 
a s for any s. Since we cannot choose sq = (Ao = oo) as in the D = 1 case, to use Eq. (|170| ) 
we must specify a renormalization prescription that allows us to fix the coupling at some 
finite non-zero value of so- We discuss this issue in detail later in this subsection. 
The equation for F$ \P,p) is given by 



,(2), 



_i_«_^)^, 2 , (0 ,o)-^ (p ,y). 



(172) 



Integrating from sq to s we obtain 



fP(p, p ') 



ds 



1 _ e -2s' p 2 p' : 
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M, 



Insisting that Fg (p,p r ) = when p,p' = we obtain 

1 



FP(p,p>) 
At third-order we have 



2 i2\ 



7 + ln(2s p 2 p 2 ) - Ei(-2s p 2 p' 2 ) 



(173) 
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2vr ds (2vr) 
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(2vr)2 as ds * lP ' Pj+ (2^)3 as ds 
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/J 2 W)+af/J 3 W), 



(175) 



where 



/fW) 
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dfe fc (2k 2 - P 2 -p' 2 ) e -2,[ P V 2 +fc 4 -fc 2 (P 2 + P ' 2 )] 



F^foAO+ifW) 



(176) 
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In the limit p,p' — » we obtain: 

1 da 



-^=a*lP (0,0) -all® (0,0), 



where 



4 3) (o,o) 



(2vr) 3 7 



dk 2k 3 e- 2sk4 



FP(0,k)+F^(k,0) 



(177) 



(178) 



Since Fg (k, 0) = Fg (0, fc) = 0, the term proportional to a 3 in Eq. ( |177 ) is zero. 
The equation for Fg (p,p) is given by 



1 dFP(p,p') 1 



(2tt) 3 ds 



- 1^(0,0) fP(p, p >) + iW(p,p'). 



TV 



(179) 



To obtain Fg (p,p) the integrals in k and s must be evaluated numerically. 
At fourth-order we obtain 

_J_da ± 1 2 dFp\p,p') 2 da ±F (2) ( a 

2tt & + (2vr)2^ d S + (2vr)2 as ds s [P ' P j 
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-a 2 lf\p,p') + al lf\p,p') + a 4 s I^(p,p') , 



(180) 



where 
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dk k (2k 2 -p 2 -p' 2 ) e -2s [p ^ + k^k^ +P -)] 



x ~FW(p, k) +F®(k,p') +FP(p,k)FP(k,p>] 
In the limit p,p' — ► we obtain: 

^^ = ^4 2 )(0,0)-a 4 /i 4 )(0,0), 



where 



/( 4 )(0,0) = — r dkk2k 2 e- 2sk * 
2tt Jo 

For dimensional reasons Eq. (182) takes the form 



Fj 3 ) ((),&)+ if) (M) 



da s B 
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(181) 
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r(4). 



where -B2 = §= and -B4 can be obtained by evaluating 7$ (0,0) (numerically) for s = 1. In 
terms of the cutoff A we obtain: 



"ST = a* * 2 a * " a* * 4 a * • 
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(185) 



Integration of Eq. ( |182| ) leads to a transcendental equation which is solved numerically in 
order to obtain the running coupling a Sj 4. 

Qualitatively, the errors are expected to be a combination of inverse powers of A and 
powers (or inverse powers) of ln(A) coming from the perturbative expansion in powers of ot\ 
for the coefficients of the irrelevant operators and the perturbative approximation for a\. 

As pointed out above, to completely determine the renormalized hamiltonian at a given 
order we need to specify the coupling at some scale Ao- The simplest way is to choose a 
value for the 'exact' a\ . Formally, this fixes the underlying theory; i.e., if we had the exact 
hamiltonian (to all orders) we could obtain the exact values for all observables. However, since 
the hamiltonian is derived perturbatively and must be truncated at some order in practical 
calculations, we can only obtain approximate cutoff-dependent results for the observables. 
Moreover, in this case the errors cannot be directly evaluated, since the exact values remain 
unknown. As an example, we calculate the bound-state energy choosing a\ = 1.45 at 
Ao = 100. In Fig. 8 we show the binding-energy as a function of the cutoff A using the 
following approximations for the interaction: 



(a) marginal operator with coupling (a\ ), 

i 2 /2\2 

(b) marginal operator with running coupling renormalized to second-order {a\^) 



V\(p,p') 



a A , 2 _(p 2 -p' 2 ) 2 



2tt 



(186) 
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(c) marginal operator plus second-order irrelevant operator with running coupling renormal- 

(2) 

ized to second-order (a\ t 2,F^ ), 



ax ' 2 + (^Y ifwo 



(d) marginal operator with running coupling renormalized to fourth-order (0^,4) 

, a A , 4 -^pi 

V\(p,p) = — — e a* ; 



(188) 



(189) 



(e) marginal operator plus second-order irrelevant operator with running coupling renormal- 

(2) 
ized to fourth-order a\^,F^ ), 

V\(p,p') = 



«A,4 / «A,4 
2vr V 2vr 



,(2) 



Fr(p, P ') 



(190) 



We see that as the approximation is improved the cutoff dependence is reduced. As 
A — > 00 all curves should approach the same binding-energy, which corresponds to the exact 
value, and as A becomes small the perturbative approximation breaks down. 

A similar prescription is to find a\ at Ao that produces a given binding-energy, Eq. Since 
the fitting is implemented using a truncated hamiltonian, this a\ is an approximation that 
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Figure 8: The binding energy for the two-dimensional delta-function potential with various ap- 
proximations for the SRG hamiltonian. The exact theory is fixed by choosing a = 1.45 at 
A = 100. 
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Figure 9: The SRG errors in the binding energy for the two-dimensional delta-function potential 
using various approximations for the SRG hamiltonian. The exact theory is fixed by choosing 
E = 1 at A = 100. 



becomes more accurate if we use a larger Ao and/or include higher order operators. Although 
in this case we can evaluate the errors, the scaling analysis becomes complicated as A — ► Ao 
because at this point we force the energy to be the exact value and so the error is zero. 

As an example, we calculate the bound-state energy when the coupling is fixed at Ao = 100 
to give Eq = 1. In Fig. 9 we show the 'errors' in the binding energy using the same 
approximations listed above for the potential. As expected, all error lines drop abruptly to 
zero when A — > Ao, where the running coupling is chosen to fit what we define to be the exact 
binding energy. Away from this point we can analyze the errors. With the hamiltonian (a) 
(unrenormalized) we obtain a dominant error that scales like In (Ao/A), corresponding to the 
leading order error. With the renormalized hamiltonian (b) the dominant errors scale like 
[In (Ao/A)] -2 , indicating the elimination of the leading order logarithmic errors. With the 
hamiltonian (c) there is a small shift, but no significant change in the error scaling. The 
added irrelevant operator may remove errors of order (Ao/A) 2 , but these are smaller than the 
remaining [In (Ao/A)] -2 errors. With hamiltonians (d) and (e) in the range of intermediate 
cutoffs (Eq << A 2 << Aq) there is also only a shift in the errors. The dips in (d) and (e) 
correspond to values of A where the binding energy equals the exact value. 

This behavior is a perturbative artifact that can be understood in the following way. 
Consider the Schrodinger equation with potential (d). Rescaling the momenta p — ► Xp we 
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As shown in Fig. 10, the coupling renormalized to fourth-order, a\4, approximately freezes 
for small A and as a consequence the bound-state energy scales like E\ ~ A 2 x constant 
eventually becoming equal to the exact value and then deviating again. With the hamiltonian 
(e) the behavior is similar, with the dip occurring at a different value of A because of the 
irrelevant operator. As in the D = 1 case, for small values of A the lines converge, indicating 
the breakdown of the perturbative expansion. 

An alternative prescription is to use the potential derived in subsection 5.1 as the starting 
point for the similarity transformation. We introduce a large momentum cutoff A, define 



Oi so =o = «A 



4tt 
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and set all irrelevant operators to zero at so = 0. Note that the coupling oa is fixed at Ao = 
oo by fitting the exact binding energy. With this definition the similarity hamiltonian with 



39 



no similarity cutoff becomes well-defined and we can set all of the similarity transformation 
boundary conditions at s = 0. In some sense, the problem becomes similar to the D = 1 
case. The previous derivation remains essentially the same. The only modification is that 
all integrals over momentum are cut off (p < A). 
At second-order we have 
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The resulting second-order running coupling and irrelevant operator are given respectively 
by 
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At third-order we have 
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In this case, 
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Since i^ A (0, 0) / the term proportional to aj? A in Eq. ( |177| ) does not vanish. To obtain 
a s ,A,3 we evaluate / s A (0,0) and solve Eq. ( |177|) numerically. 
The equation for F^(p,p') is given by 
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(3) 
To obtain F^(p,p) the integrals over /c and s must be evaluated numerically. In the limit 

sA 4 — > 00 with s fixed at some non-zero value 
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and a Sj A becomes indeterminate, requiring the coupling to be fixed at some s$ 7^ 0. In this 
way, we recover the result of the previous prescription. 

Although less trivial, this prescription allows a more transparent error analysis. We can 
also extend the calculation to larger values of the similarity cutoff, A, and analyze the errors 
in a different scaling regime. In Fig. 11 we show the errors in the binding-energy obtained 
using the following approximations for the potential with A = 50: 



(a) marginal operator with coupling (a\ ) 
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(b) marginal operator with running coupling renormalized to second-order (a\ a 2)1 
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(c) marginal operator plus second-order irrelevant operator with running coupling renormal- 
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Figure 11: The SRG errors in the binding energy for the two-dimensional delta- function potential 
using various approximations for the similarity hamiltonian. The exact theory is fixed by regu- 
lating the "bare hamiltonian" using a sharp momentum cutoff, A, and letting the bare coupling 
depend on A such that the binding energy is fixed. We use A = 50 and E Q = 1. 
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There are clearly two distinct scaling regions when an additional large momentum cutoff 
A is placed on the initial matrix and a similarity cutoff is then applied. When the similarity 
cutoff is larger than A, we see scaling behavior that is similar to EFT. Curves (a), (b), and 
(d) all have the same slope. None of these hamiltonians contains irrelevant operators, but 
the marginal coupling differs in each. All results become exact as the similarity cutoff goes 
to infinity, and these curves are close to one another because the coupling runs little in this 
region. Curves (c) and (e) show that there is a power-law improvement when irrelevant 
operators are added, and that once again when the similarity cutoff is larger than A an 
improvement in the running coupling makes little difference. These results closely resemble 
those in EFT and we conclude that it is the leading irrelevant operator that leads to this 
improvement. Even though the coupling in front of this operator is approximated by the 
first term in an expansion in powers of the running coupling, the coupling is sufficiently 
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Figure 12: The SRG running coupling for the two-dimensional delta-function potential renormal- 
ized to second and third-order obtained with aA =oo — ► «a = 47r/ln ( 1 + ^- 



small that this approximation works well and the operator eliminates most of the leading 
power-law error in curves (a), (b), and (d). 

When the similarity cutoff become smaller than A we see a crossover to a more complicated 
scaling regime that resembles the SRG scaling discussed above. The error displayed by curve 
(a) approaches 100%, while the running coupling introduced in curve (b) reduces the error 
to an inverse logarithm. Improving the running coupling in curve (d) further reduces the 
error, and we see that curve (c) crosses curve (d) at a point where improving the running 
coupling becomes more important than adding irrelevant operators. As above, the best 
results require us to both improve the running coupling by adding third-order corrections and 
add the second-order irrelevant operators. In no case do we achieve power-law improvement, 
because as we have discussed there are always residual inverse logarithmic errors. Had we 
fit the running coupling to data, as we would do in a realistic calculation, we would obtain 
power-law improvement and the residual error would be proportional to an inverse power of 
the cutoff times an inverse power of the logarithm of the cutoff. 

In Fig. 12 we show the running coupling at 2nd and 3rd order. Although the 3rd order 
corrections are small for all A and vanish when A — ► oo, the improvement resulting from this 
correction in Fig. 11 is significant. 

The scaling behavior with a large momentum cutoff A in place is complicated, but it 
is fairly straightforward to understand it and to find a sequence of approximations that 
systematically improve the non-perturbative results. The calculations become increasingly 
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Figure 13: The SRG running coupling for the two-dimensional delta-function potential renormal- 
ized to second-order obtained with a\ 0=O0 — > «a = 47r/ln ( 1 + -|- 



complicated, but at each order one must improve the running coupling, or fit it to data, and 
add higher order irrelevant operators. In a field theory we need to let A — ► oo and study 
the scaling behavior of the theory in the regime where A<A. Although we do not display 
a compete set of figures, in Fig. 13 we show what happens to the running coupling as A is 
increased, with the bound state energy fixed at one. 

As is evident in the exact solution, as A increases the coupling decreases. When A>A, 
the coupling runs slowly and stays near its asymptotic value. As A approaches A the coupling 
begins to run noticeably, and when A becomes much less than A the coupling approaches a 
universal curve that is insensitive to its asymptotic value. Plots of the error in the binding 
energy for various approximations and different values of A closely resemble Fig. 11, with 
two scaling regimes whose boundary is A = A. 

We close this section by reminding the reader that in all of these calculations there is only 
one free parameter. In a realistic calculation we would fit this parameter to a binding energy 
and we would expect to see residual errors in other observables that is inversely proportional 
to powers of the cutoff and logarithms of the cutoff. 
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6 Conclusions 

We have illustrated the EFT and SRG methods for producing effective cutoff hamiltonians 
using the one-dimensional and two-dimensional delta- function potentials. We have shown 
that the addition of irrelevant operators with couplings tuned to fit low energy scattering 
observables leads to errors that scale as inverse powers of the cutoff in EFT. We have shown 
that the SRG with coupling coherence leads to errors that scale as inverse powers of the 
cutoff and logarithms of the cutoff. EFT typically leads to smaller errors than the SRG, 
but the cost is an increased number of parameters that must be fit to data. The SRG with 
coupling coherence requires the same number of parameters as the underlying 'fundamental' 
theory, but the cost is exponentially increasing algebraic complexity to remove errors that 
contain inverse powers of logarithms of the cutoff. 

Our examples do not adequately illustrate the real utility of these renormalization tools, 
because our examples are easily solved with other methods. EFT can be used when the 
underlying theory cannot be solved and even when it is not known. All that is required is 
a knowledge of the degrees of freedom and symmetries that appear at 'low' energies. The 
SRG can also be employed when the underlying theory is not known, but in general it offers 
no advantages over EFT unless that effective theory is valid over very many energy scales. 
When the effective theory is the same as the fundamental theory, but with cutoffs introduced, 
or when it applies over very many scales, it may be possible to fix some of the couplings in 
terms of the fundamental couplings without using extra data. 
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The leading-order and next-to-leading-order EFT expansion parameters 
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